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VLBI detection of an AGN pair in the binary black hole 
candidate SDSS J1536+0441 

M. Bondi^ and M-A. Perez- Torres^ 
ABSTRACT 

We present first pc-scale radio imaging of the radio-quiet candidate binary 
black hole system SDSS J1536+0441. The observations were carried out by 
the European VLBI Network at the frequency of 5 GHz and allowed to image 
SDSS J1536+0441 with a resolution of ~ 10 mas (~ 50 pc). Two compact 
radio cores are detected at the position of the kpc-scale components VLA-A and 
VLA-B, proving the presence of two compact active nuclei with radio luminosity 
Ln ~ 10'^° erg s~^, thus ruling out the possibility that the two radio sources are 
both powered by one 0.1 pc binary black hole. From a comparison with published 
8.5 GHz fiux densities we derived an estimate of the radio spectral index of the 
two pc-scale core. Both cores have flat or inverted spectral index, and at least for 
the case of VLA-A we can rule out the possibility that synchrotron self-absorption 
is responsible for the inverted radio spectrum. We suggest that thermal free-free 
emission from an X-ray heated disk wind may be powering the radio emission in 
VLA-A. 

Subject headings: quasars: individual (SDSS J153636. 22+044127.0) — galaxies: 
active — radio continuum: general 



1. Introduction 

SDSS J153636.22+044127.0 (hereafter SDSS J1536+0441) is certain ly a peculiar and 



compl ex object. The attention to this low redshift QSO was drawn by iBoroson fc Lauer 



(I2OO9I . hereafter BL09) in a search for quasars exhibiting components with multiple redshifts 
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among the Sloan Digital Sky Survey (lAbazajian et al.l 120091 ) . They found three line systems 
at different redshifts: two sets of broad-line emission at z = 0.3889 and z = 0.3727, and 
a third set of narrow absorption lines at z = 0.3878. Narrow emission lines were found 
associated only to the higher redshift system. BL09 interpreted these peculiar features in 
terms of a massive black hole binary (BHB) system within the same galaxy with separation 
of ~ 0.1 pc and masses of Mq^ and Mg^, each with its own broad line region and sharing 
the same narrow line region. BHB systems are expected as the results of galaxy mergers, but 
few compe lling candidates have survived scrutiny. Furthermo re, the so-called "final parsec 
problem" (IBegelman et al.l Il980l : iMilosavljevic &: MerritI l200ll ) is not currently understood. 
In particular, dynamical friction with the stellar background is ineffective in shrinking the 
binary below separations smaller than 1 pc where gravitational radiation can complete the 
coalescence of the two black holes within a Hubble time. 

After the BL09 publication, SDSS J1536+044 1 was target of s e veral observations with 
the aim to confirm or dispute the BHB scenario. IWrobel fc Laoii (120091 . hereafter WL09) 



used the VLA to image the quasar at 8.5 GHz and detected two radio emitting components 
separated by 0.97" (5.1 kpcj^. This suggests the possibility that SDSS J1536+0441 is not 
a BHB system separated by ~ 0.1 pc but a double quasar separated by ~ 5 kpc and 
probably residing i n a moderately rich cl uster of galaxies. High r esolution ground based and 
HST observations ( iDecarli et al.ll2009al Jbl: iLauer &: Borosonll2009l ) indeed detected an optical 
counterpart to the secondary radio component but were not conclusive on the nature of the 
optical counterpart itself (obscured AGN or elliptical galaxy) and on its relation with the 
peculiar optical spectru m. The scenario of a supe rposition of two AGNs was dismissed by 
Lauer fc BorosonI (120091 ) and IChornock et al.l (120101 ) on the basis that the regions responsible 
for the two sets of broad lines are spatially coincident, even if some doubts remain given 
that this result is based on observations ca rried out with a see ing larger than the angular 
separation of the two objects along the slit ( IDecarli et al.ll2009bl ). 



A third possible explanation was raised by lChornock et al.l ( 120091 . l2010l ). Based on Palo- 



mar and Keck optical spectra they conclude that SDSS J1536- j-0441 is an unusual raember 
' """" ' ^ 1^ ^ - ' cers (e.g. iHalpern &: Filippenkd Il988 



Eracleous & Halpern 


1994 




Strateva et al. 


2003 




pendentlv bv ICaskell 


( 


2010 


). The main problem 



striking differences betw een the Balmer-line pro file in SDS S J1536-I-0441 an d all t he other 



double-peaked emitters (ILauer fc BorosonI 120091 ). Finally, iTang fc Grindlayl ( l2009l ) argued 



that SDSS J1536+0441 could be both a double-peaked emitter and a BHB system. 



^We assume a flat cosmology with Hq = 71 km s ^ Mpc ^ and ilm — 0.27. 
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All these pos sible scenarios, together with the ejected black hole hypothesis, were criti- 
cally discussed in iLauer fc Borosoru (120091 ) . No matter what is the real explanation for the 
observed properties of SDSS J1536+0441, what we can surely say is that this is a puzzling 
and possibly unique object. 

In this Letter, we present for the first time high-resolution (~ 10 mas corresponding 
to ~ 50 pc) images of the 5 GHz continuum radio emission of SDSS J1536+0441, obtained 
with the European VLSI Network (EVN). These data detect pc-scale radio cores both in 
SDSS J1536-I-0441 and in the companion at 0.97" east of the quasar proving the presence of 
two AGNs in this region. 



2. Pc-scale Imaging 

The European VLB! Network (EVN) was used to observe SDSS J1536+0441 at 5 GHz 
on 2009 October 23. The observations were carried at 1024 Mbit s~^ sustained bit rate to 
exploit the large bandwidth capabilities of the EVN, with an array which included all the 
EVN antennas. The pointing position was centered on the VLA-A component detected by 
WL09. Observations were phase-referenced to the calibrator J1539+0430 with a duty cicle 
of 5 minutes. The total time on source for SDSS J1536+0441 was about 5 hours. The strong 
and compact sources J1613+3412 and J0154+4743 were used as fringe finders and bandpass 
calibrators. 

The data were correlated and calibrated at the JIVE correlator. Standard a priori 
gain calibration was performed using the measured gains and system temperatures of each 
antenna. The amplitude calibration was refined using the phase reference source. Further 
data inspection, flagging and imaging were performed by the authors using the NRAO AIPS 
software. No polarization or self-calibrations were performed. 

The la r.m.s. noise in the final images is about 15 /xJy. 



3. Results and Discussion 

Figure 1 shows the field of SDSS J1536-I-0441 imaged with the EVN at 5 GHz using 
natural weighting, which resulted in a beam of 12 x 7 mas in position angle 10° and a la r.m.s 
noise of 15 /xJy beam^^. Both radio sources, VLA-A and VLA-B following the notation used 
by WL09, are clearly detected with signal-to-noise ratios (SNR) of 50 and 15, respectively. 
Contour plots of VLA-A and VLA-B are shown in Fig. 2. 
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We also imaged the two radio sources with a shghtly better resolution (7x6 mas) 
but worse sensitivity, using a different u-v data weighting function. These images are not 
shown here but both set of images were used to find consistent fitted Gaussian parameters 
for VLA-A and VLA-B. The sources were fitted with two-dimensional elliptical Gaussians. 
The flux densities, positions and errors are: for VLA-A, 5* = 0.72 ± 0.06 mJy, a( J2000) = 
15^36'^36^2232, 5(72000) = +04°41'27".069, and ctvlbi = 0.003 mas; for VLA-B, S = 0.24 ± 
0.03 mJy, a(J2000) = 15^36'"36^2881, 5(J2000) = +04°41'27".054, and ctvlbi = 0.003 mas. 
VLA-A appears slightly resolved when fitted with a single elliptical Gaussian component. 
The deconvolved fitted sizes are 6'vla-a — 3.2x2.5 mas with an estimated error of 1 mas. The 
high SNR for this component provides a high level of confidence for the deconvolved fitted 
size. There is some indication that also VLA-B could be slightly resolved with deconvolved 
size of about 2 mas, but in this case given the limited SNR this value should be considered 
as an upper hmit. 

The positions are in excellent agreement with those determined by WL09. The separa- 
tion between the two components is confirmed to be 0.97". 

One of the goals of the proposed EVN observations was to determine the origin of the 
radio emission in VLA-B and in particular if it could be associated with VLA-A, e.g. as 
being a compact jet/hot-spot or a mini-lobe ejected by VLA-A. No extended emission is 
detected in the region between the two compact radio sources. Given the compactness of 
VLA-B on the pc-scale we can now confirm without any doubt that VLA-A and VLA-B are 
two compact AGNs and we can rule out the possibility that the two radio sources are both 
powered by a 0.1 pc binary (see WL09). 

Comparing our flux densities at 5 GHz with those obtained by WL09 we can derive 
an estimate of the radio spectral index between these two frequencies. It certainly can be 
speculative to draw conclusions on the radio spectral index of VLA-A and VLA-B based 
on observations made at different resolutions and epochs but it is a flrst step, waiting for 
more appropriate observations, and it is illustrative of the general trend of the spectral 
properties. As far as radio variability is concerned, the VLA 8.5 GHz and the VLBI 5 GHz 
observations are separated by 8 months and significant variability (e.g. ~ 20% with a few 
cases wi th larger variations) is observed on such time range in radio-quiet and radio-loud 



quasars (IBarvainis et al.l 120051 ). Keeping in mind this limitation, it is worth noting that the 
fiux density of VLA-A at 5 GHz, 0.72 mJy, is lower than the 1.17 mJy measured at 8.5 GHz, 
while the fiux density of VLA-B is rather constant: 0.27 mJy at 8.5 GHz compared with 
0.24 mJy at 5 GHz. Both sources would have fiat or inverted radio spectra, and VLA-A 
in particular would exhibit a strongly rising radio spectrum. WL09 already suggested that 
VLA-A could have an inverted radio spectrum to explain the non detection of the radio 
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source in the 1.4 GHz FIRST sky survey (jWhite et al.l 119971 ) and our 5 GHz observations 
confirm this possibihty. Another possibihty is that we are missing flux in the 5 GHz VLBI 
observations. If there is some extended emission (on the scale of ~ 100 mas) our VLBI 
observations might not be deep enough or have the adequate u-v coverage to detect it. This 
possibility is rather improbable since a significant amount of extended, and therefore steep 
spectrum emission at 5 GHz, should have been detected in the 1.4 GHz First survey given the 
1 mJy threshold of the VLA survey. Therefore, in the remaining discussion we will assume 
that any missing flux is not significantly affecting the derived spectral indices. 

With the measured fiux densities at 5 and 8.5 GHz, VLA-A would have a spectral index 
a ~ —0.9 (with S^u) oc z/~" ). Such an inverted spectrum, even if found in a few radio quie t 
quasars, is quite unusual ( iBarvainis et al.l Il996t iKukula et al.l Il998t lUlvestad et al.l l2005l ). 
Synchrotron self-absorption is the most invoked cause to explain the fiat or inverted spectra 
in radio-loud or radio-quiet AGNs. For self-absorption to occur, the brightness temperature 
must be comparable to the ki netic temperature of t he synchrotron electrons. The brightness 
temperature in Kelvin is (e.g. lUlvestad et al.ll2005l ): 
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which for VLA-A gives = 9 x 10*^ K and for VLA-B Tf, > 6 x 10^ K (assuming un 
upper limit for the size of VLA-B of 2 x 2 mas) . These brightness temperatures, at least for 
VLA-A for which we have a measured fi tted size, are too low to affect the spectra unless the 
magnetic fields are unrealistically high (iGallimore et al.lll996l ). 



The optical counterpart of VLA-B is very red (jPecarli et al.l l2009bl : iLauer fc Boroson 



20091 ) and with the detection of a pc-scale fiat spectrum radio core with observed radio 



luminosity = vLy at 5 GHz of 0.6 x 10 erg s~ is best interpreted as an obscured AGN 
rather than an elliptical galaxy. The brightness temperature of VLA-B is only a lower limit 
since the size is not constrained and we cannot rule out the possibility of synchrotron self- 
absorption affecting the flat radio spectrum as well as the contribution of thermal free-free 
absorption/emission. 

The low brightness temperature of VLA-A rules out synchrotron self-absorption as the 
origin of the inverted radio spectrum. Free-free absorption should be discarded as well, since 
the view to the AGN is unobscured in VLA-A where we see optical continuum emission 
from the AGN and the broad-line region. We suggest that the radio emission from VLA-A 



could be interpreted as thermal free-free emission from a disk wind (IGallimore et al.l Il996 



Blundell &: Kuncid 120071 ). The disk wind is heated by th e X-ray continuum an d therefore we 
can expect a link between the radio and X-ray emission. iPanessa et al.l (|2007| ) found a linear 
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correlation between radio {Lr at 5 GHz) and X-ray {Lx) core luminosities in an optically 
selected sa mple of nearby Seyfe rt. This result was confirmed and extended to radio-quiet 
quasars by iLaor fc Beharl (120081 ) using the carefully selected and almost complete Palomar- 
Green quasar sample. They found that L^jLx ~ 10~^ where is the radio luminosity 
at 5 GHz and Lx is the bolometric 0.2-20 keV X-ray luminosity. Q uite remarkably, the 
same correlation, known as Giidel-Benz relation (jGiidel fc Beng|l993l ). holds for coronally 
active stars, therefore covering a range of abou t 15 orders of magnitude in luminosity. WL09 
used the X-ray luminosity measured by Swijt (lArzoumananian et al.ll2009l ) and extrapolated 
the 5 GHz luminosity from their 8.5 GHz measurement assuming a radio spectral index 
a = 0.5, obtaining L^/Lx = 5.9 x 10~^. Having a direct measure of the core flux density 
at 5 GHz we can refine this measurement. The observed monochromatic luminosity at 5 
GHz of VLA-A is Lsghz = 3.8 x 10^° erg s~^ Hz~^, and the radio luminosity at 5 GHz is 
Lji = uL^QHz = 1-9 X 10^° erg s~^ which gives L^/Lx = 1.4 x 10~^, assuming all the X-ray 
emission is associated to VLA-A. 



These VLBI observations were not meant to solve the so-far unsatisfactory interpretation 
of the puzzling properties of SDSS J1536-I-0441. The debate is still open about all the possible 
scenarios summarised in the Introduction. What is now clear is that both VLA-A and VLA- 
B are powered by their own AGN, and this should be taken into account in future analysis. 
Further radio observations are necessary to confirm the spectral shape of the radio cores and 
therefore the origin of the radio emission. 



4. Summary 

We have presented the first VLBI pc-scale imaging of the candidate binary black hole 
system SDSS J1536-I-0441 observed by the European VLBI Network at 5 GHz. Both the 
VLA-A component, associated with SDSS J1536+0441, and the companion object, VLA-B, 
at 0.97" east of the quasar are detected with high SNR (50 and 15 respectively). The two 
radio nuclei appears barely resolved and no extended larger scale emission is detected. The 
main results can be summarised as follow: 

• We detect a flat spectrum pc-scale radio nucleus at the position of VLA-B confirming 
that both VLA-A and VLA-B are powered by their own AGN. Given the radio and 
optical properties, VLA-B is most likely associated to an obscured AGN rather than a 
passive elliptical galaxy. 

• At the position of VLA-A we detect a slightly resolved radio nucleus with a strongly 
rising spectrum with a ^ —0.9 between 5 and 8.5 GHz. We rule out synchroron self- 
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absorption as the cause of the inverted radio spectrum because the derived brightness 
temperature of VLA-A is too low. We suggest that thermal free-free emission from a 
disk wind provides the simplest explanation for the inverted radio spectrum in VLA-A. 

We derive a value of Ln/Lx = 1.4 x 10~^ for VLA-A that is to tally consisten t with 
correlation found f or radio-quiet quasars and Seyfert galaxies (jPanessa et al.l 120071 : 
Laor fc Behailboosh . 
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Fig. 1. — EVN grey-scale image at 5 GHz of the SDSS 1536+0441 field. This image was 
obtained applying natural weighting and has a resolution of 12 x 7 mas at position angle 
10° VLA-A (on the right) and VLA-B are clearly distinguishable with SNR of 50 and 15 
respectively. 
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Fig. 2.— Contour plots of component VLA-A (right) and VLA-B (left) of SDSS 1536+0441. 
The beam size is 12 x 7 mas at position angle 10°. The la r.m.s noise is about 0.015 mJy 
beam~^. First contour is 3 times the noise and each inner contour is -\/2 times brighter. 



